Universality in Dynamical Formation of Entanglement for Quantum Chaos 
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Dynamical formation of entanglement is studied for quantum chaotic bi-particle systems. We 
find that statistical properties of the Schmidt eigenvalues for strong chaos are well described by 
the random matrix theory of the Laguerre ensemble. This implies that entanglement formation for 
quantum chaos has universal properties, and does not depend on specific aspects of the systems. 
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At the dawn of the quantum mechanics, entangled 
states of spatially separated bi-particle systems triggered 
fierce discussions because of their paradoxical property 
PJ. These discussions lead to the finding of the inequal- 
ities which measure how quantum systems deviate from 
the classical reality |2j. Recently, a renewed interest has 
been paid to entanglement since it is recognized as an 
invaluable resource in quantum computation and quan- 
tum communication 3]. Thus, the information theory 
and the circuit theory for entangled qubits has become 
fast growing fields. 

On the other hand, in the quantum theory of open 
systems, dynamical formation of entanglement between 
a system and its environment is a target of the research. 
In these studies, it is supposed that uncontrolled forma- 
tion of entanglement causes the loss of quantum coher- 
ence in the system and leads to the appearance of the 
classical reality. Then, the loss of quantum coherence 
is considered to provide the foundation of the statistical 
physics Also in the quantum computation, dynam- 

ical formation of entanglement between qubits and their 
environment is a crucial subject since it can give rise to 
erroneous operations of quantum computers. 

In the traditional study on the dynamical formation 
of entanglement, the system is supposed to interact with 
an environment composed of a large number of degrees 
of freedom. Then, the interaction with the environment 
causes complicated uncontrollable formation of entangle- 
ment, leading to the loss of coherence in the system. 

However, interaction with an environment is not the 
only cause of complicated processes of forming entangle- 
ment. Actually, dynamical formation of entanglement 
by quantum chaos has been an active topic these days 



000 In their studies, quantum chaos gives rise to 
increase of the von Neumann entropy (or a modified ver- 
sion of it) to the extent that the state approaches the 
maximally entangled state, which is an equally weighted 
superposition of entangled bases. Therefore, quantum 
chaotic systems can also be the origin of the loss of co- 
herence, even when they are of a small number of degrees 
of freedom. 

This indicates that there are several causes where com- 
plicated processes of forming entanglement takes place. 
Then, the next question is how we can classify, if pos- 
sible, those processes of forming complicated entangle- 
ment. Are there any universality classes according to 
which we can differentiate seemingly similar processes of 
forming entanglement ? 

In this letter, we show that there actually exists a uni- 
versal class in the dynamical processes of forming en- 
tangled states. Here, we investigate the time evolution 
of Schmidt eigenvalues for coupled systems of quantum 
chaos, and indicate that, for strongly chaotic systems, the 
distribution of Schmidt eigenvalues approaches a class 
of random matrices, i.e., the Laguerre unitary ensemble 
0- So far, the Laguerre unitary ensemble has received 
only a little attention in the theory of random matrices 
[ljj,[n|. We will discuss possible implications of this uni- 
versal class, and suggest its extensions for other cases. 

We consider a wave function of a bi-particle system, 
<$>(xi,X2)- For simplicity, we assume that the dimensions 
of the Hilbert spaces for each degree of freedom are the 
same, and will be denoted by N. Then, the wave function 
is represented as 
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Here{|« >i \i = 1,2, . . . , N} and{|j > 2 \j = l,2,...,N} 
are the orthogonal bases for each degree of freedom, and 
Aij(i,j = 1,2,..., AT) are complex- valued coefficients. 
For the N x N matrix A whose components are Ay, we 
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perform the singular value decomposition [T^. 

A = C/Al/" 1 (2) 

where f and V are N x N unitary matrices, and A is a 
diagonal one whose components are Aj(i = 1,2, ... , N). 
Here, we can always make \ non-negative real and set 
Xi > Xj for i < j. We transform the bases from {\i >i} 
and {\j > 2 } to {\i >i} and {\j > 2 } using U and V, 
respectively so that the wave function is described as the 
Schmidt decomposition form, 



N 



$(xi,x 2 ) = ^A i |?>i (g>|i> 5 



(3) 



For the representation J3J , the von Neumann entropy for 
the reduced density matrix, in which one of the degree of 
freedom, x\ or x 2 is traced out, is easily evaluated as 



N 



S = -^Afln(A|). 



(4) 



We note that only the diagonal matrix A contributes to 
S. 

From now on, we construct a random matrix theory 
for the coefficient matrices, A of wave functions. First, 
we assume that the occurrence probability for the ma- 
trix A is a function of the independent components A^. 
Moreover, we assume that the probability is invariant 
under local unitary transformations. These assumptions 
determine the joint probability function P({A»j}) as 



P{{A ij })^dA ij = CeM-^M^A)]]JdA 

i,3 i,3 



(5) 



where a and C are constants. Transforming the argu- 
ments of P from Aij into Xi and the parameters which 
determine U and V , and integrating the probability with 
respect to the variables other than A, , we obtain the prob- 
ability function with respect to the eigenvalues 



PL({s i })l[de t ^C'Ulsi-Ej 



JJdei, (6) 



i<j 



where we introduced the new variables Si(= N 2 Xf), 
which take the values from to oo, and the normaliza- 
tion constant C @E3- F rom the joint probability iJBJ, 
we can get the n-point correlation function as 



y* oo />oo 

/ •••/ P L (e 1 ,...,e N )de n+1 ---de N 

JQ JO 



(N - n)\ 
det[K N (ei, ej)]i t j-i t ... tn , 



where 



N-l 



K N (ei,£j) = ^ <Pk(£i)tpk(£j) 



(7) 



(8) 



Here, we take ipk{£) = e^ £ / 2 Ck{e) with £fe(e) being the 
A;-th normalized Laguerre polynomial, hence we call the 
distribution © the Laguerre unitary ensemble. From 
Eq. J7|), we obtain the one- level density 0] and the two- 
level cluster function [111 as 



<r N (e) = Ri(e) = K N (e,e) 



(9) 



and 



T 2 (e 1 ,e 2 ) = -R 2 (ei,£ 2 ) + Ri(£i)Ri(£2) - [K N (e lt e 2 )] 2 , 

(10) 

respectively. In order to discuss the statistical property 
of the level spacing, we need to unfold the level spectrum 
[To|. Thus we introduce the rescaled level ui{e) as follows 



L0(S) 



<jN(e')de'. 



(11) 



According to Nagao and Slevin |l5j], we also introduce 
the renormalized two-level cluster function as 



T 2 (uj,uj') 



r 2 ( £ M, £ (a/)) 

Vi?i(eH)iii(e(w'))' 



(12) 



In comparing the Laguerre unitary ensemble with an 
ensemble of wavef unctions, we will divide the distribu- 
tion of the Schmidt eigenvalues into the following three 
regions: The hard edge, i.e., those eigenvalues lying near 
zero, the soft edge, i.e., those eigenvalues lying near the 
largest one, and the bulk region, i.e., those lying away 
from the hard and soft edges. Note that specific features 
of the Laguerre unitary ensemble emerge only in the re- 
gion where eigenvalues are near zero, i.e., the hard edge. 
For the other regions, it is known that characteristics of 
the Laguerre unitary ensemble are similar to those of the 
Gaussian unitary ensemble |l5|. 

In the following, we study the dynamical formation 
of entanglement for a coupled system of two degrees of 
freedom, and compare the Schmidt eigenvalues of wave- 
functions to the Laguerre unitary ensemble . We give the 
Hamiltonian as 

H T = Hi(x\,pi) + H 2 (x 2 ,p 2 ) + Hi 2 (pi,p 2 ), (13) 

where 



Hi(xi,pi) = 2sin 2 (^) +kiJ2 s Mxi)S{t-n) (i = 1,2), 

n=l 

(14) 

Hi 2 (pi,p 2 ) = 4c pp sin(y)sin(y). (15) 

Here fej (i = 1, 2) are the parameters for the kick strength 
and c pp is a coupling constant. In order to make the 
dimension of the Hilbert space finite N, we impose the 
periodic boundary conditions not only for Xi but also 
for pi, i.e., < Xi < 2ir and — it < pi < ir (i = 1,2). 
These conditions require that the Planck constant be h = 
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FIG. 1: Time evolution of S and A; . The time evolution of 
S, Af, A| and A§ is plotted for the strong chaos with TV = 
128. The initial state is a product of the coherent states 
centered at (xi,pi) = (ir/2,n/4) (i = 1,2). The increment 
of S is slowly saturated before reaching the maximal entropy 
S max = ln(JV). 



2n/N, and we choose N = 128. For comparison between 
strong chaos and weak chaos, we choose the parameters 
as follows. For strong chaos, their values are ki = 3.0, 
&2 = 2.5 and, for weak chaos, they are k\ — 0.7, k 2 = 0.2. 
In both cases, we choose c pp = 0.05. 

As seen in Fig. 2] the time evolution of Q by the 
strong chaos shows a temporal generation of entangle- 
ment between x\ and x%, where we choose as an ini- 
tial state a product of two coherent states. We notice 
that the variation of the entropy consists of two stages. 
While the entropy increases rapidly in the initial stage, it 
starts to saturate for the latter stage. In the initial stage, 
the largest Schmidt eigenvalue decreases rapidly with the 
other ones growing from zero. For the latter stage, these 
Schmidt eigenvalues start to avoid each other Jjjl| . The 
same behavior was found in a similar system 14]. Such 
avoided crossings remind us of the behavior of the en- 
ergy eigenvalues for quantum chaotic systems. For weak 
chaos, saturation in the increase of entropy also takes 
place. 

Both for weak chaos and strong chaos, we construct 
an ensemble of 10 6 wavefunctions, respectively, by col- 
lecting them after the entropy starts to saturate. From 
the ensembles of the numerical data, we evaluate the one- 
level density R\ (e) and the two- level correlation function 
R2(s,e') directly, and estimate the renormalized two- 
level cluster function Tb(u;, a/) through Eq. 112|) . 

Adachi and Iida 0] compared the one-level density of 
Schmidt eigenvalues in the coupled standard map with 
the prediction of the random matrix theory, and found 
that they agree well with each other. We confirmed their 
results for our system [13j . In order to see higher order 
statistics, we compare, in Fig. [21 the renormalized two- 
level cluster functions T^w, a/) between the Laguerre 
unitary ensemble and the ensembles produced by the 
time evolution. For the hard edge, T 2 (0,w) is depicted 
(Fig.JSJa)). For the strong chaos, the renormalized two- 



level cluster function agrees with that of the random ma- 
trix theory, although the fluctuations are not negligible. 
For the weak chaos, it behaves similarly to that of the 
random matrix theory, while the deviation from the ran- 
dom matrix theory is remarkable inw< 1. For the bulk 
region, in Fig. H^b), T 2 (oj' , oj' + uj) is averaged over cu' in 
the bulk region. For the strong chaos, the cluster function 
agrees well with that of the random matrix theory. How- 
ever, the cluster function for the weak chaos takes nega- 
tive values after it falls. This indicates that the Schmidt 
eigenvalues show weaker avoided crossings. For the soft 
edge, f 2 (N - w, N) is depicted (Fig. I^c)). For both the 
strong and weak chaos, the renormalized two- level cluster 
functions do not agree with that of the random matrix 
theory. In particular, an oscillation with large amplitude 
appears for the weak chaos. 

In Fig. the distribution of level spacing for the bulk 
region is shown for the Gaussian unitary ensemble, the 
ensemble of the strong chaos, and that of the weak chaos. 
Here, we can see that the ensemble of the strong chaos 
agrees well with the Gaussian unitary ensemble. On the 
other hand, for the weak chaos, the level spacing distri- 
bution shifts to smaller level spacings than for the strong 
chaos. This is consistent with the results of the renormal- 
ized cluster function showing weaker avoided crossings. 

We have shown that, when the coupled system exhibits 
strong quantum chaos, the set of the Schmidt eigenval- 
ues approaches the Laguerre unitary ensemble. When 
the coupled system shows weak quantum chaos, the set 
deviates from the the Laguerre unitary ensemble. This 
indicates that the Schmidt eigenvalues for coupled sys- 
tems exhibiting quantum chaos belong to a universality 
class of the Laguerre unitary ensemble. 

The universality of the Schmidt eigenvalues implies 
that their distribution will be insensitive to how the cou- 
pled system is divided into two subsystems interacting 
with each other. In other words, the boundary between 
subsystems has the flexibility that the properties of sub- 
systems do not depend significantly on how we define the 
boundary. This is exactly what we expect when the cou- 
pled system behaves in a classically statistical way. Thus, 
the universality of the Schmidt eigenvalues can be con- 
sidered as one of the crucial features of the effects called 
decoherence. 

However, the approach to the the Laguerre unitary 
ensemble is not complete even for the coupled system 
exhibiting strong quantum chaos. This is shown by 
the slight deviation around the soft edge. This means 
that dynamics within the Hilbert space spanned by these 
Schmidt eigenvectors still behave in a nonstatistical man- 
ner. Then, the following questions are the next target of 
our study: How does the boundary between those be- 
having statistically and those behaving in a nonstatisti- 
cal way change in time ? Is there a natural boundary 
where we can separate the coupled system into the two 
subsystems ? How does the boundary differ depending 
on the choice of initial conditions ? 

There exist a couple of ways to extend the idea of uni- 
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FIG. 3: Level spacing distribution. The distributions for the 
level spacing Sl in the bulk region are plotted. 

versality classes into other cases. One possibility is to 
consider the coupled system where the symmetry exists 
between the two subsystems. We expect that the ex- 



istence of the symmetry will lead to other universality 
classes. Another way of extending our idea is to consider 
the coupled system where the sizes of the Hilbert spaces 
for the subsystems differ significantly |14j . This situation 
corresponds to the cases where the traditional studies of 
the statistical physics and the measurement theory con- 
sider. One possibility is that the universality classes can 
be described by the extended Laguerre ensemble. If so, 
we can think of the following problem: How does the de- 
coherence effects vary between the cases when the system 
interacts with the other one of the same size and those 
when it interacts with the environment? 

The study into the above questions is under progress 
and the results will be published elsewhere. 
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